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We have studied primordial non-Gaussian features through bispectrum and trispectrum analysis 
from a model of potential driven DBI Galileon inflation originating from background supergravity 
and Gauss-Bonnet terms. We have explicitly shown the violation of the widely accepted Mal- 
dacena theorem and standard Suyama-Yamaguchi relation in squeezed limit configuration which 
leads to the resolution of the well-known sensitivity problem between the non-Gaussian parameters 
(/ivi, Tjvi, pjvl) and tensor to scalar ratio (r). Our analysis thus overcomes a generic drawback of the 
wide class of DBI inflationary models which was, of late, facing tension from observational ground. 
Hence large primordial non-Gaussianities can be obtained even from single field DBI Galileon and 
■ hence these class of models can be directly confronted with the forthcoming results of PLANCK. 
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I. INTRODUCTION 



The physics of the early universe is a very rich area of theoretical physics, for there is a plethora of potential models 
that solve, at least partially, the well-known problems of the standard cosmological paradigm. Inflationary cosmology 
is the most successful branch which addressed all of these problems meticulously. This can however be explained by 
several class of models originated from a proper field theoretic or particle physics framework. But from observational 
point view a big issue may crop up in model discrimination and also in the removal of the degeneracy of cosmological 
parameters obtained from Cosmic Microwave Background (CMB) observations [|], §]. In this context the study of 
primordial non-Gaussian feature acts as a powerful computational tool to discriminate and to rule out several class 
of proposed inflationary models. In the very recent days the analysis of bispectrum and trispectrum derived from the 
study of primordial features of non-Gaussianity from different models of inflation has thus become an intriguing 

aspect in the context of inflationary model building as well as studies of GMB physics. 
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The DBI inflationary model [|9| , one of the most interesting possibilities to realize large non-Gaussianity of 

the cosmic microwave background (CMB) temperature fluctuations. But this class of models has a serious drawback, 
commonly known as generic sensitivity problem which states that the non-Gaussian parameters derived from bispec- 
trum and trispectrum analysis is highly sensitive to the tensor to scalar ratio . This results in a tension between 

these two parameters from observational ground. To be specific, the parameter space of non-Gaussianity and tensor 
to scalar ratio as obtained from DBI models contradicts the future predictions of PLANCK [|| which has already 
started giving early results - the science run results are soon to follow. Very recently, a natural extension to the wide 
class of DBI models has been brought forward by tagging Galileon with the good old DBI model, resulting in so-called 
"DBI Galileon" p5|-p7| which, amidst its successes, results in unwanted degrees of freedom like ghosts, Laplacian and 
Tachyonic instabilities. In a recent paper the present authors have demonstrated how these unwanted debris can 
be removed in spite of keeping all the good features of |2^ intact by proposing DBI Galileon in D3 brane in the 
background of ^=1,15=4 SUGRA derived from D4 brane in 7V=2,I?=5 bulk SUGRA background. Starting from this 
model. In the present paper, our prime objective is to carry forward this rich structure of DBI Galileon as proposed 
in in order to resolve the above mentioned sensitivity problem by studying bispectrum and trispectrum calculated 
from three and four point correlation function originated through cosmological perturbation, working upon third and 
fourth order perturbative action in De Sitter (DS) and Beyond De Sitter (BDS) limit. Subsequently, we demonstrate 
that it is indeed possible to have a parameter space for both non-Gaussianity and tensor-to-scalar ratio (r) consistent 
with observational predictions of PLANCK. 

The plan of the paper is as follows. First we exhaustively study primordial non-Gaussian features from the third 
order perturbative action through the non- linear parameter /^vl calculated from bispectrum (in non-local and equi- 
lateral limit configuration) including all possible scalar - tensor type of interactions in the different polarizing modes. 
Hence from the fourth order perturbative action we derive the expression for other two non-linear parameters g^L and 
tnl through trispectrum analysis. Finally, we have explicitly shown that the violation of the well known Maldacena 
theorem [ p9) , pO| and standard Suyama-Yamaguchi relation pl[ , in the squeezed limit configuration leading to 
the solution of the sensitivity problem of non-Gaussian parameters. 

II. THE BACKGROUND MODEL AND PERTURBATIVE ACTION 

For systematic development of the formalism, let us briefiy review from our previous paper [^8| how one can 
construct the effective 4D inflationary potential for DBI Galileon starting from M ~ 2,2? = 5 SUGRA along with 
Gauss-Bonnet correction in the bulk geometry and D4 brane setup leads to an effective N ~ 1,1) — 4 SUGRA in the 
D3 brane. Here the total five dimensional model is described by the following action 



where 



.(5) _ "(5) 
■"5 



c(5) _ 0(5) I c(5) I 0(5) I c(5) (0 T\ 

^Total — '^EH ~^ ^GB "l" DA brane BulkSugra V^'"^^ 



Sfn = ^ / d'x^/^) [i?(5) - 2A5] , (2.2) 
I d^.v^ - 4i?-^^(^)4') + rU . (2.3) 

ZK5 J L J 



The D4 brane can be decomposed into two parts as 



^D4 brane ~ '-'DBI ' '-'WZ' 

where the DBI action and the Wess-Zumino action is given by 



0(5) _ c(5) I c,(5) /r, 

•-'D4 brane ~ "^DBI ^ "^WZ' \'^-^) 



SdIi = d^^eM-"^)]/- (7(5) -|-27ra'^^(5)), (2.5) 



^WZ = - f / E„=0,2,4 '^n A exp (^2 -|- 2na' U form 



T!a CijFcD 
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(2.6) 



§^^BijBKLdG^'^dD<^^ + ^^BijFcd) + 2n^a^TiCoFABFcD - Tiv{^)\ } 
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where T(4) is the D4 brane tension, a is theRegge Slope determines length scale of the string, exp(— $) is the closed 
string dilaton and Cq is the Axion. In eqn(2.1) M — 2,7) — b bulk supergravity action can be written as 



o(5) 
Bulk Sugra 



1*- ~ r™"9V-'If» 



2 I 2 " " • 'I- ^ q 

Fermionic + Chern ~ Simons + Pauli mass] 



yc.f!{D,nr)iD'''r) (2.7) 



In the present context 5-dimensional coordinates = {x" , y) , where y parameterizes the extra dimension com- 
pactified on the closed interval [—ttR, +TrR] and Z2 orbifolding symmetry is imposed. 

Applying dimensional reduction technique the total effective model for D3 DBI Galileon in background M=\, 2?=4 
SUGRA is described by the following action: 



(4) 



+ h {Cil)R 



«/375(4)^(4) 



a/^(4) n(4) 



al3 



A{6)R^ 



(4) 



+ ^3 



(2., 



where 



K{c^,X) 
h 



p 



1 - 2QXf - Qi 



C5G{^,X) - 2XAf(T,rt) - 



1 + ^ (241(2) - 24^(9) - 16^(10)) 



"(4)C(2) 



7^ 



(2.9) 



All the constants appearing in equation(^ 
Coleman Weinberg potential is given by: 

1/(0) = 



and equation(|2.9|) are explicitly mentioned in [p8l . The one-loop corrected 



9^ (24C(24) - 144Z(4) - 64^(5) + 144^(7) + 64^(8) + 192^(11)) - 2.p4> ^(l 



1 



i2m 



(2.10) 



, C2 = (^l3RI{2)f2-gv^Z{T,T^)A{l3)), C4 



— '-^ — — ) and Dq — 0. Here the energy scale of inflation has a GUT scale window 0.658 x 



where Cq = (T.ii>a + l3RI{2)fo + Z{T,T^)A{l3)v'^y C_4 = T^va 

(liRT{2)h 

lO^^Gey < ^ < 0.667 x W^GeV. 

Using this model our next goal is to exhaustively study the primordial features of non-Gaussianities by deriving the 
explicit form of bispectrum and trispectrum in DS and EDS limit which can solve the generic problem of well known 
DBI inflation in our framework. Moreover the explicit calculation confirms that effective DBI Galileon in D3 brane 
including the Gauss-Bonnet correction is necessary ingredient to resolve the sensitivity problem of DBI inflation. To 
serve this purpose we begin with our discussion to the following perturbative action which can be expressed as: 

/ 



/ 



dt d^x 



JC-CC + jC-hh + J-lvv 



Chh 



Cfhhh ~\~ ^hhv ~\~ ^hvv H~ C^vvv ~l~ -^CC^ ~^ ^C,vv ~\~ ^C^hv 



\j}ower spectrum contribution 



bispectrum contribution 



-c. 



cccc 



+ ^CCCh + ^CChh + ^ 



-Chhh 



"Qhhv 



trispectrum contribution 

(2.11) 

where Cab-^ ^abc and Cabcd represent second, third and fourth order Lagrangian density respectively with (a, 6, c, d := 
C,v,h). In equation( ^.ll ) (, v and h stands for the primordial curvature (scalar), vector and tensor perturbation 
modes. It is imperative to mention that, in this context, we neglect all the contributions from vector modes in all 
order of perturbation theory. Further, we neglect the contribution of tensor modes in the fourth order perturbative 
action as well as in the four point correlation functions. To explore the non-Gaussian features from our model we 
start with the following generalized ansatz for the primordial curvature perturbation and tensor perturbation: 

A'"(f) X^ix) + ^ (/^TO"'™ {(Xcm^ iXGix)\ - ((A'G(f)), (Xcixm + ^ (g'^TT'^"' {Xcix))^ {XciS)), (^'^(f)) 

+ ^ (/i^r')"' {{Xam, {Xg{x)), {Xg{x))^ {Xg{x)), - ((A'g(x)), {Xg{S)\ (^'^(f)), {Xg{x)),}} 

+ - {X{x)) 

(2.12) 
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where {a,b,c,d = C,h). Specifically (A'gCS'))^ = J -j^j (^Xcik)^ ^exp (j,k.x^ where 
hij{k)elj^\k). Moreover in this context the n-point correlation functions are given by, 

( :n=l,3,5 V(a,6,c,d) 

{{XG{x))^{Xc{x))f^{XG{x))^{XGix))^ n terms) = < finite{^0) :n=4,6,8....'^{a,b,c,d) and n=2 with a^b 

I :n=2 with a=b. 

(2.13) 



III. TREE LEVEL BISPECTRUM ANALYSIS 
A. Three scalar correlation 



To calculate the scalar bispectrum for D3 DBI Galileon we consider here the third order perturbative action up to 
total derivatives. Using the uniform field gauge analysis the third order perturbative action for three scalar interaction 
can be written as: 



where 



j-{a'YsC)~aYscld\ 



can be calculated from the second order action |28 



'CC 



dtd^xa^Ys 



Here C^{i = 1,2,3, 
Ys 



are dimensionless co-efficients defined as: 



MpL 



2t 



3 %r- 3 ' ^ 



HYs 



Ys Ys 

Li{Liai + as) + ai2 + (ag + Lia4)— + — 

ti ci 



3M2^ 



l + ^l(aL,{Ys-t,}) 
a at 

,C4 = -^0 + 2ti 



1 ^PL 


d 
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dt 





d (A. 



dt \ tj 



t^ 
''1 



t^ 

''1 



(3.1) 
(3.2) 

(3.3) 



(3.4) 



6 2ti & dt^ ' " " 1^ ti 

and the co-efficient of involves spatial and time derivative in equation(5.9) is defined by the following expression: 



7^ = ^ {{d.OidkX) - d~^d,d, mQ){d,x)]] +PiCC - ^ {(50' - d-^d,d, [(5.C)(a,C)]} . (3.5) 
In this context ^3 — 2I5, — —^iMem. and 7?. — > as fc — > at large scale. Throughout our analysis we use the 



Hubble parameter as if = 



A(4)-H87rG(4)V(0) 
91 



from our previous paper 



Additionally Li 



Mi 



\HMj,j^~4,XGx , 

—picj, and X — '(^sO- Now following the prescription of in-in formalism in the interacting picture the three point 
correlation Junction both for DS and BDS, after some trivial algebra took: 



aki)ak2)ah)AH[%) 



8 „0 

= {2nf5^^\ki + k2 + k3)Bccc{ki,k2, h), 



ccc 



|0) 



(3.6) 
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where the total Hamiltonian in the interaction picture can be expressed in terms of the third order Lagrangian density 
as (i?i„t(77))^^^ = = - / d^x (jOs)^^^. Moreover the bispectrum B(^,^(^{ki,k2,k3) is defined as: 



(27r)47'2 



6 



(3.7) 



where the symbol ; 1 is used for three scalar correlation. Here ^^^^(fci, k2, k^) is the shape function for bispectrum 
and is used for normalization of E-mode polarization expressed in terms of the new combination of the cyclic 
permutations of two-point correlation functions given by 



Pi = Pdkr)Pi{k2) + Pdk2)Pdk3) + Pdk3)Pdki)- 



(3.8) 



In this context /]^£"' represents the non linear parameter carrying the signature of primordial non-Gaussianities of 
the curvature perturbation in bispectrum. The explicit form of characterizing the bispectrum can be expressed 

as: 



r f local] 
IJNL;1\ ds 



Ci / 2 



2 3 1 
+ X! ^i^i ~ Ir^ X! 



3 1 ^ 



kikj 



4c2 I 2 



3 -, 3 

ij'(«>j) = l i,j(i=Aj) = i 

3 ^ ^ 
~ 9 X! ^i^i ~ kik2k3 ^ fcifcj 



3i?C'3 {kik2k3f ^ YsC4 



2MpL 

3 



8M|,^ 



« J'(«>j) = l 



3 



1=1 

^ 3CeH^ (fcifc2/c3)' ^ Cr^f' 



1 , , , 3fcifc2A;3 



+ 



^clMl^K^ 



2 

2 ^ -' 

^fcifc2fc3E^^-^fci^2fc3i^'-6 E fcf^l-Efcf + ^E^' 

i=l = l 1=1 1=1 



3 



(3.9) 



r Aocal~\ 
L/jVZ,;lJ 



10 



BDS 



kik2k3 



3(l-ey-s^) 



— 1 



r(i^.) 



r(i) 



Ci 



+ 



+ 



3(1 - ey - s 
Ys 



2Ys 



S\2 



J^sMnc - 1) + -ils{i>) 



J^eleinc - 1) + ^26(z>) 



+ ^^4(2^) + ^^2^5(1^) 

,S\2 ' 



(3.10) 



C,(l-ev.-4)2 C3(l-ev.-5^ 

+ 7:77-^ 2:7(1/) + —= 18( 



8c2 



where 




Ii{x) 


= Cos ( [x — i 


Mx) 


= Cos ( [x — i 


Ia{x) 


= Cos ( [x — 5 


Xc,{x) 


= Cos ( [a; — 5 



2j 2)r(i + ^) 



2+XY^ pt.2_l+KY^ 7,21,3 J/ 



X] = " Cos ([x ■ 



1] 7r^^ (6+x)r(3+x) 
12 ' 



T^-^E.>,fc.fc.-^fclfc2fc3 ,I,(x) = i^ 



(ktk^k^f r(3+x) 



Cos([x-i] f), 



It{x) = Cos ( [x - 
Ts{x) = Cos ([x 



+ i^i^ [(3 + .x)r(i + x) - r(2 + x)f ] } 
^[r(i + x) + r(2 + x)^] + 

[r(i + x) + r(2 + x)^]} 



11 1) 2±2 [r(i + x) + r(2 + x) (fcifci+i^+Mfci + (3 + x)^i|^)] |Iii:fM + I^k^ + {^^-\^ ^ 



2J 2/ 



j^] Tl\ j (fci.fc2)fc3 



(3 + x)r(l + x) + (3 + x)r(2 + x)f - r(3 + x)-^ + 1^1^ [(3 + x)r(l + x) 



+ (3 + x)r(2 + x)f - r(3 + x)|i] + i^k^ j(3 + ^)r(i + ^) + (3 + ^)r(2 + x)^ - r(3 + x)^] } 



(3.11) 
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Here we have defined K = ki + k2 + k^, x = (ri(^ — 1, i/) and 



2ev 



^3 



1 — ev ~ s 
Ys{l + Ys) 



1 + ev 



^ Ys-ev + {l + Ys)p3 ' 
1 + ey 



H 



£3 ■■= 



= I P3 

71 



Pi 

n2 



(3.12) 



Ysjl + Ys) 
1 + ev 
2{1 + Ys)^ 

2p4{l+Ys) 



2T4 - 

Ys-ey 
1 + Ys 



1 + Ys 
1 + ev 

+ P3 



(1 - 2pi) 



with four new constants pa, P4, Ta, Tl- For the numerical estimation we have further used the equilateral configuration 
(ki — k2 — kj, = k and K — ik) in which the non- linear parameter fiqL can be simplified to the following form as: 



pequil 
JnL;1 



DS 



85 
324 



1 



1 

lOe, 



1 



12 



4^T/'(0)M 



85 _ 
54' 

PL 



305 
162 



Sgx 



27c? 27c? 



85 
108c? 



{(s + Vs 



24) 



81Sg 
65 



5e 



/20Ac 



65 



V 81e. 



162c2es 



72e, 



5e 



108c? 



lOt^G 
81 



(3.13) 



requil 
JNL-l 



10 



1 



—1 



BDS 9k^ V 54 



^ir''(r.c-l)- 



r( 



sov 



y2 



2YsS^y 



3-ey /^l 
4c? 



l + ev 



S ^2 



3(1 -ey- 4) 



2Ys 



^313 



"«™'(nc-l) + %l3"«"'(:>) 



j-equil,~ 



3(1 - ey - 4) 
Is 



4c2 V 4(1 + ey) 



J-6X|«"'(nc-l) + ^I|'""(i>) 



6 q-equil / ~\ 



(1 - ey - 4)2(1 + ys)2(y5 - ey)^,,„,; , . 



(1 + Ys){Ys - ey)(l - ey - 4)^egmi , ~ 



2Ysc1{l + evY 



4c2(l + ey)2 

In the equilateral configuration for BDS 



(^) 



^equil, 



"XT' 



■equil 



(x) 



2{l-x) ' 
r^«"'(a;) = f 1 + |) Xg"""' 



(3.14) 
(3.15) 



Now using the tensor-to-scalar ratio 



r{k)] 



DS 



16e,,c., 



(3.16) 



ik)] 



BDS 



1 - ey - 4 ^ ^ 
1 — ey — s 



(3.17) 



the sound speed c^ can be eliminated from the equation( 3.13) and equation(3.14) also. Most significantly the numerical 
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FIG. 1: |[ajVariation of the /jvl vs sound speed for scalar modes eg, 'b) Variation of the /jvl vs tensor to scalar ratio r in DS 
and BDS hmit. 



value of f^^^l in the equilateral limit is obtained from our set up as 2.8 < f^^^l < 7 in DS limit and 4 < Z^^™' ^ 
in EDS' limit within the window for tensor-to-scalar ratio 0.213 < r < 0.250 |Q. This is extremely interesting result 
as it is different from other class of DBI models. The most impressing fact is that the upper bounds of /^^™' in both 
the limits are within the future observational bound as predicted by PLANCK |Q . The graphical behavior of /^^™' 
with respect to the sound speed (cs) and tensor to scalar ratio (r) are plotted in the figure (|l|) (a) and figure (|^)(b). It 
is evident from the figure(|l|) that /^^™' is not much sensitive to the tensor to scalar ratio (r) or the sound speed (cg) 
within a specified range applicable for our model in DS and BDS limit analysis. 

B. One scalar two tensor correlation 

The one scalar and two tensor interaction for D3 DBI Galileon can be represented in uniform gauge by the following 
third order perturbative action: 



Chh 



dtcPxSQhh 



(3.18) 



where the dimensionful coefficients = 1, 2 7) are defined as: 

HLiYt Yt d /Li 



1 - 



3 dt \(3j. 



T2 = Yscl, T, 



h^^iY:^^Kxx]+'2a 



Ys ^ HLiYt Y.^ 



HY, 



Tr, = 2aYTLi 



1 



^6 



AaY^ 
' Yt 



-2Y,. 



2Y4 — 



d ( aL\ 



dt V^tCt/ ' 



(3.19) 



where we use a = (pXG^x- The extra term in the action represented in equation(3.18) can be expressed as 

S(hh 



(3.20) 
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with the following gauge fixing condition: 



0. 



(3.21) 

(3.22) 



The contribution from the last term in the above action can be gauged away by making use of the following gauge 
transformation equations along with the gauge fixing condition for tensor and scalar modes: 



LiYt 



(3.23) 



Now following the prescription of in-in formalism in the interaction picture three point one scalar two tensor correlation 
function both for DS and BDS can be expressed in the following form: 



q=l-'-°° L ^ 

= (2^)3^(3) ^ ^ {Bchh}ij,u (^1' ^2, fca), 



ij;kl 



Chh 



|o) 



(3.24) 



where the total Hamiltonian in the interaction picture can be expressed in terms of the third order Lagrangian 



density as ([H^ntiv)],,.^),,, = E,'=i ( \hI%)] . = - J d^x \{£s)^^^ 



ij;kl 



. Moreover the cross bispectrum 



{Bchhjij-ki {ki,k2,ks) is defined as: 



lli=l "^i 



(3.25) 



where the symbol ; 2 is stands for one scalar two tensor correlation. Here {^Chh)ij.i^i (^1,^2) ^s) is the shape function 
for bispectrum and the polarization indices are u = 1{E — mode),2{E^B — mode),3{B — mode). We adopt the 
following normalization depending on the polarization in which we are interested in: 



Pdki)Pcik2) + Pcik2)Pc{k3) + Pi{k3)Pi{ki) :u^l(E-mode) 
Pc(ki)Ph{k2) + Pdk2)Ph{k3) + Pc{kz)Ph{ki) :u=2(E®B mode) 
Ph{ki)Phik2) + Ph{k2)Ph{kz) + Ph{ki)Ph{ki) :u=3(B-mode). 



(3.26) 



Consequently [//^£°2] represents the non-linear parameter which carries the signature of primordial non- 
Gaussianities of the one scalar two tensor interaction. The explicit form of [/^£"2] "ij-ki characterizing the bispectrum 
can be calculated as: 



ij',kl 
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7S,(fci,fc3,fc2)]"a(^'Hfci,fc3,fc2)} 



(3.27) 



where in the DS limit the characteristic polarized coefficients in the tensor basis are given by: 



iij,kl 



'ij,kl 



,mn 
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'ij,kl 



- u 
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Hj^kl 



^ ^ klmkln 



= kl 



' '^3^kl 



(3.28) 
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with polarization index u = 1{E),2{E 1^ B),3{B) and Mij-ki = (Ea ei^j(^)4;^^(^))- 

In this context the momentum 

dependent functions can be written as: 



1 c\.klkl{cski + K) 



1 clkl + 2clcTkl{k2 + ks) + 2c^c^ki{k^ + fczfcg + fc|) + 4(^2 + k3){k^ + k^k^ + kl) 



1 cl(^klkl{K + CTks) 



K 

g{4) _ 2 clc^k\k2k1 ^(5) _ 2d^k 2k^{'iCsk\ + K) 



^(6) = g(4)^ g(7)=^(5) 

where we define K_ := Cgfci + CT{k2 + fcs). In BDS limit we have 



(3.29) 



r flocall " 
[JNL;2\ ij. 



kl 



^ I 22i/s+4l/t-12 



^(^ 
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(3.30) 



where in the BDS Umit 

^Jpiki,k2,k3) 
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(^4)"j.fc; 
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ij;kl 



k'i^ik^k^YT 
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(I-i/t) c 

(fci.fe) yi23 + (fci.fe) 



k^^^ikik^y^ 

4 



+ 



Jp{k3,k2,ki) 



ij;kl 



r(7 + p-4z/T-2i^,,; 



fc3"(fc2fcl)''^ 

r(9 + 7; - 4//:;' - 2/^,) 



4i/T- 



2y, 



+ ( Ara.fc^ ) I213 + fe.A^il I231 + (ki-h) I312 + ffc; 



■'^ij,kh 



Iki) 



Y:- 



321 



ij,kl 



J123 + 1/132 + 1/213 + >/231 + 1/312 + 1/321 



A/"" 



C123 + Cl32 + C213 + C231 + C312 + C321 



ij^kl 



W123 + Wl32 + W213 + W23I + W312 + W32I 



A/"" 



[klmki^' {X123 + X132} + k2mk2^' {X23I + X213} + fcsrnfta^' {X312 + X321}] Klkl^ 



mn.m n 



(3.31) 



with 



i/i (^o) kb, kc) 

1/3(^0) ^6) ^c) 
Siikai kbi kc) 

1/5(^0) kb,kc) 



= -^ij kl il - >^t) , J2{ka, h, kc) = Jl {K, kb, kc) KCs ' ^ 

il,r{,L L _i Zj jKl L _i ZJ^Kl 

= Klki (i - ^t) [(fc^ + kl + kakb) + (I - i/t) kaikb + kc)] 
= Kl kl (1 - "t) [kbkcikb + kc) + ka{kl + kl + hkc)] , 

a r _ ^ ^ -| tt 

~ -^ii kl l^h^c "I" kakbkcikb "I" ^c) I ) Jsika, kb, kc) — iJ^Ij klk^kb^c^ 

ij;kl L Aij;kl ■" 

{includes 3 permutations of a, b, c), 



(3.32) 
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^abc 



T{8 + q- AvT - 2vs) 



abc 

9 2u,-§-§ AuT-^-if 



j^8+q—4l^T- 



Hi = H2 = Li, H3 



abc _ Oi^5^sC^ /3 



Hi {kah + kbkc + kcka) 



Ha 



Cg Crp 



{with a,6,c= 1,2, 3 with a ^ b ^ c) , 



ti V 2 



v. 

2 



,Af' = Cs ''crj.'^KAf'^Af" 



kakbkc [^-Vs] +kl{kb + fee) 



) AT", 



A abc 



a^YsC^ 2 
7 kg^k^kc, 



(3.33) 



T{Q+p-AvT-2vs) 



-Avt — '^^s 



af' = (I - i^r) - [ka{kb + fee) + kl + kl + kbkc] + Q - j^t) A:^, 



- - Z^,s I , 



) "2 



a6c n^^^f^ 3 



(Z-l^ Cg C'jn -/il^ . 



r.abc 



^abc 



klikb + fcc) ( ^ - j + {fca(A;6 + A:^ + fcbfcc) + kbkeikb + kc)} - i^Tj - 



I ~ vt\ kl(kl + kl + fcfefce) + ( ^ - 



k\k\ 



\ ^ ( ^ ' kakbkc{kb + kc) 



ka kbkc 



^ - I^T j A;b/Cc + f ^ - I's ) fca(^6 + kc) 



abc _ _uIU^U^ 
) "7 ^ '^a'^b'^ct 



(3.34) 



^ [jp(fca,fcfa,V,J r(7+p-4z.T-2:/.) r 7^ « f ^ 

= > , ,, , , — — z 5- — ta {includes 6 permuations of a,o,c) 

^r,kl j^^ kl'{kbk,y^ ^2,..-f-1^4.r-t-f ^7+p-4.r-2., ^ ^ 3 , , ) 



V5; 1 / ,2 r(8 - AvT ~ 2v,) , g^nc^ ^ r(6 + p - 4:.t - 2i/,) 



kl'{kbkcY^ \ "cf=-t4'^^-'#^8-4-T-2., ^ P ^2.,-§-2^4.r-^-4^6+p-4.r-2.3 

r(7 + p-4z/T-2z/«) 



kl'(kbkcY-^ cf^-S-ic^^T-f -^^7+p-4^^-2^3 ' 



(3.35) 

For both DS and iJfS' limit the overall normalization factor for three types of polarization can be expressed as: 

{8 :u=l(E-mode) 
128 :u=2(E(^B mode) (3.36) 
2048 :u=3(B-mode). 

Further, to make the computation simpler without loosing any essential information we reduce the complete set in 
terms of the two-polarization (helicity) mode instead of four complicated tensor indices. For this purpose let us define 
a reduced physical quantity: 

0'(fc) = /i,,(fc)elf (3.37) 
in terms of which the one scalar two tensor correlation is defined as: 

(C(fci)0^'(fc2)0^'(fc3)) = {2T:)H{k+k + k)B\^i:t{\kuk2M- (3-38) 
where the cross reduced hispectrum is defined as: 

D(^2;A3)/jr r r ^ _ {^irtVl {X^-\^) _ 6 rj-ioca/n^K^^lAs) p2 /„ ^a-. 

lli=l '^i 
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Here applying the basis transformation the explicit form of [/]^£"2] ^^^'^"^ characterizing the crossed bispectrum can 
be calculated as: 



r f local 
[TnL;2 



2^2 



(3.40) 



where the DS polarized coefficients in the helicity basis ( 712^3(^1) ^^2, ^3) ) can be evaluated as 



(3) 



1 6 A/g 

ki ■ ks 



[kl-{X2k2 + X3k3f]\ 



kl 



^ A2 ,A3 



' A2,A3 



A2 jAs 



k2 ■ k3 



' A2,A3 



(4) 



^(1) 
' A2 jAs 



(5) 
^A2 jAs 



' A2 jAs 



= kl 



A2 jAs 



K 



(fci - - fc3)(fci + A;2 - A;3)(fci - + ^3) [fci - (A2fc2 + A3A;3)^]" , 



^A2 jAs 



and in BDS limit 



(3.41) 



[JnL;2\ 



Cs C 



BDS 3j:i,kf 

^ I 221/3 +4i^T- 12 



r(^^.) 



r(f) 



r(j^T) 



r(f) 



X 



(fci)"- {k2k3r 
(l-ev-4)^(l~ey^.^)Vi(^) ^ 



32.Fi (Vi)"'^^'^=' + 4J-2 (V2)"'^^'^^ + 2 (V3)'" 
+ h {Vir^''^' +h (V5)"'^==^^ +Te (Ve)"'^^'^^ + (Vr! 



A2 ;A3 



where in the BDS case 
[jp{k 



(Vi) 



m;A2;A3 



E 

p=i 



^;A2;A3 



fcr=(A;2fc3)''^ 



■ + 



Jp{k2,ki,k3) 



k^^ikik^r^ 



+ 



Jp{k3,k2,ki) 



>;A2;A3 " 



(V2) 



m;A2;A3 



r, \ (fcl.fe2) I (fc2-fc3) I (fc3-fcl) ^ 

_ [fei=(fc2fe3)''T k^'iklkayT "I" feJ'*(feife2)''Tj r(9+p-4z/T - 2z/s) 



kl'{k2kiY^ 



(I - H 



A^i.A^;) ri23 + [ki.k^ Yrs2 + [h-h) Y213 + [h-h^ Y231 + [h.h) Y312 + [h-ki) 



21/3 -f -3 4i/T- 



(V3) 



m;A2;A3 



— ~^ j^9+p—iuT — 2us 

Y321 



(3.42) 



^i;A2;A3 



(V4) 



«;A2;A3 ^ 2 i I - ly. 



(Vs)"'^''^' 
(Ve)"'^''^' ^ 2 

|-^^^«;A2;A3 



J123 + i7l32 + i7213 + J23I + J312 + J32I 

m;A2;A3 



C123 + C132 + C213 + C231 + C312 + C321 

VV123 + Wl32 + W213 + W231 + W312 + W32I 



^A2A3 



7«;A2;A3 f y i V \ _L 7";A2;A3 f y _i_ Y \ _l 7"|A2;A3 f y i V \ 
^1 \^123 + ^132 / + ^2 \^231 + ^213/ + ^3 \^312+^321/ 



(3.43) 
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with 



Jsika, kb, kc) 

yli(kai kh, kc) 
yJbika J kb, kc) 



u:\2;\3 



u;\2 ;A3 



u;\2 :A3 



2(i-^T) 



yJlikai kh, kc 



u:\2-.\3 



yJl (ka, kh, kc) 



2A2A3 (I - vt) [{kl + kl + kakh) + (I - Vt) ka{kh + kc)] 



= 2(1 



i/t) [X^khkcikh + kc) + \lka{kl + kl + khkc)] 



(3.44) 



= 2X1X1 [klkl + kahkcikh + kc)] , Jaika, h, kc) 



— 2z Ao Ao /ct^ 



(includes 3 permutations of a,b,c), 



and 



a 6c 



E 



j^j7J,(fca, fcb, fee) 



'";A2 ;A3 



r(7 + p- 4z^T - 2z^s 



_ 2p _ 14 



ka" {khkc)""^ ^2i/,-f-I^4iyr-^-^^7+p-4<.r-2,., 



In the equilateral limit we have 



u;(A2;A3) 



[includes 6 permuations of a,b,c), 



(3.45) 



r f local 
UNL;2icqml 



DS 



y2 2 



+ ^0 



jil\^{k,k,k)Ygip)ik,k,k)} 



(3.46) 



flocall uK-^2;A3) 



r flocall ^ 
L/JVL;2J , 



BDS 



lOQ^" (i - ^t)' (jcs + 2c^)fc)^-^+^-^-^ [Co. ([.. - 1] f )] ^ [Co. ([.T - k] f )] • 



r(i) 



r(:/T) 



r(i; 



^ (l-ev-4)^(l-ev-.|:)Vi(^) ' 



+2(j-3(V3):_ 

\'u;A2;A3 , -J^ /v- 7 \'i'';A2;A3 



where each coefficients and functions appearing in equation((^ll|)-(|05|)) are evaluated in equilateral limit 



(3.47) 



C. Two scalar one tensor correlation 



The interactions involving one tensor and two scalars are given by the following third order perturbative action: 

+ 3^6i9^/i»jV',*V'j} + / dtd^xScch 



/■ 3 3 f 3^1 3^2 • • 3^4 2 3^5 2 

= dtd X a < —hijC,iC,j + —hijC,tC,] + yshtjCiiJjj + —a hijC,tip,j + —o hjC.^Cj 



(3.48) 



where the last term of the equation(3.48) can be expressed as: 



LiKxx 



Licr 



(3.49) 



LiYs (kxx o 
Yt \ 2 + ' - 



00 



8^- 
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and the dimensionful coefficients = 1,2 7) are defined as: 



3^1 = Y,cl, 

y2 



3^5 



LiKxx /^r 2 \r 2\ , T v2 
V^Cs - YtC^) + LxY4 



1 , HL.Kxx 



-2 + 



+ 



4 



+ 2HLiYTa-YT-(L,a), 



3^3 = n 

3^4 = n 



3 ^ d / KxxLi ^ a \ _ Ljj^Yt\ I KxxLi ^ a 



2 dt \ 2 Yt 

(Yr-Iyxx4)L 



2 

XL 

AYt 



(3.50) 



This last term of the perturbative action can be gauged away by making use of the following gauge transformations 
for the tensor and scalar modes: 



hii + 4 



LiYs iKxx ^ o 



L\Yt 



LiKxx 



YtLI (kxx a 



CiC: 



(3.51) 
(3.52) 



Following the prescription of in-in formalism in the interacting picture three point two scalar one tensor correlation 
function both for DS and BDS can be expressed in the following form: 



7 -.0 

= {2nfS^^\h+k2 + k3){Bcch}ki {kiMM). 



aki)ak2)hM{h)\[Hl^t{v) 



|0) 



(3.53) 



where the total Hamiltonian in the interaction picture can be expressed in terms of the third order Lagrangian density 
as {[Hint{v)]ki)cch = Eg=i il^intiv)] ) = -Jd^x [(A)^^ J ^ • Here the cross bispectrum {B^ch}ki is defined 



as: 



(3.54) 



where {Ac^c,h)^.i is the two scalar one tensor correlation shape function and the symbol ;3 represents two scalar one 
tensor correlation. Consequently the non-linear parameter [/^£'?|] can be expressed as: 



r rlocal 1 
IJnL;3\ 



kl 



ily<i { [7!?']" (fci, fe, fe)/C(«' (fci, fe, fca) + 



where the polarization coefficients in DS' limit are given by: 



(fe,^i,fc3)/C(«)(fc2,A;i,fc3)} 
(3.55) 



7iP 



kikhiJ^ij^i^s), 



(2) 

Tki 



(6)" 

Tki 



^3 
'>-l"'2 



(3) 

Tki 



ko 



1^ 



(4) 

Tki 



u 1,2 
_ ^ 



(3.56) 
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and the momentum dependent functions are given by: 

(1) 1 clih + fc2)(fcf + fcifca + kl) + 2c^CT(fcf + fcifca + k^)ks + 2c^4(fci + fca)^! + 4fc| 

t^(2) 1 4^3[2c^(A;? + 3fciA:2 + kl) + 3csCT(/ei + ^2)^3 + 4A;i] 
^ = F F ' 



1 clc\.klkl{cski +^) 

1 c^A:|[c;^(A-i + A-,)(2A:i+A:2j+ 3c,r, (2A-i + A',) A-, + 2r^.A:2] 
if K 



3 



/C(5' = [c^(fci + A;2)(fc? + 3A;iA;2 + fcf ) + ^clcrikl + 3fcifc2 + A;|)fc3 + 4c,4(fci + fcs)^! + 4fcf] , 

^(6) _ 1 C^fcffcK^ + CTfca) 

with ^ := Cs(A:i + ^2) + CtA:3. 
In the BDS hmit 



[JNL-Mkl 



_ lOL^^^Af.rM K'-^^"^^-' [Cos{[u, - 1] f )] ^ [Cos {[ut - |] f )] ' 



SOS 3E?=ifc| 

^ I 24!^s+2i/T-12 



r(i) 



4^, n2/- <^n4 



r(i) 



(l-ev-4)"(l-ev-.^)ViW ^ 
YiYT44~glMl,^ 



where in BDS limit we get: 

(v,). 



(3.58) 



{k2ik3j + k3ik2j) _^ {kiikaj + fcsjfcij) (fcHfc2j + fc2»fcij) 

ki;^{kik^Y^ k^/i^kik^Y' 



fVs) = 

(Vi) = 

\ / ij 



k\^{k2k^Y 
3 

vt 



fVe) 



k\ 



{k2tk3jPl23 + kzik2jP\32) (fcl»fc3jP2i3 + ks^kij P231) {kiik2j P312 + k2ikijP32l) 

k'i^{k2k3Y^ k^^ikiksY^ k^^{kik2Y' 

{k2ik3jRi23 + k3ik2jRl32) {kuk3jR2l3 + k3ikijR23l) {kiik2jR3\2 + k2ikijR32l) 

kr{k2k3Y' k2'{kik3Y' k^'^{kik2Y'- 

k2ik3jRi23 + k3ik2jRl32 ] ( fcljfc3j-R213 + k3ikijR231 ] ( ^14^2^-^312 + k2ik\jR32l ) 

/ _|_ ^2 \ / J 1,2 \ / 



fcr(fc2fc3)''^ 



2 {htk:ij + A;i,:A-2j) _^ ^^2 + /':itfclj) _^ ^^2 (fcl»fc2j + ^'i^fclj) 



k^-'{kik3Y^ 

■ + kl 



+ kl 



fcr(^i^2)"= 



k'i^{k2k3Y^ k2^{kik3Y' k^^{kik2Y 

( k2ik3jLi23 + k3ik2jLi32 ] ( kuk3jL213 + k3ikijL231 ] ( kuk2jL3i2 + k2ikijL32i 
V ^ _|_ iu2_\ /_ _(_ 1^2 \ 

fcr(fc2fc3)^= ' k^-^{kik3Y^ 



+ ^3- 



k^-^ikikiY' 



(.3.59) 



with 



, p=l I 



r(9+p-4i^, -2i/t) 



4i/, 



-6 2iyT- 



(3.60) 



_ 2 _ _!_ 

Oi = 1,02 = «Cs ^Cj,^K,03 = - {kah + kbkc + kcka),OA = -ikahkc 
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Pahc 



r(8+p-4!/T-2!/,) 



2i/T- 



i j^S+p-ivs -2vt 



abc r^abc 



--vt] {kah + hkc + kcka) + fc„ 



(3.61) 



^-vt] kakbkc + kl{kb + fee) 



^ abc 



ka kbkcj 



Rabc = Li(--Ut)J2-^: 



p=l 



^abc 



+ 



ti 



A;? 



yJlikai kb, kc) 

J^aika, kb, kc) 

J^iikai kbi kc) 
Sbikai kb, kc) 

and 



1 abc 



ics ^c~^KAf',Af^ 



[/V — ♦ — * — * "1 ^ 
Jq{ka, kb, kc) 

9=1 

--vt) [kah + kbkc + kcka + kl] 



r(7 + p-2j/T-4i^.) 




-^5 — kgkbkc 

J2{ka,kb,kc) 



iCs ^ Crp^ K 



^ 1 (^a J kb, kc) 



vt {kah + hkc + kcka} 



kjtc + o - kcka + 7T - k^kb 



klkl +[i:-fs] klhkc +[i:-^t] kahk. 



--VT\kckb + kakbkc 



yli(ka, kb, kc) 



Seika, h, kc) 



zkf,khk^ 



(3.62) 



Rabc = klLiO + 



Labc — L\0 



Lia^Yg (I - Ug) 



p=i 



,abc 



r{8 + p- Ai^s - 'ii'T) 



4Ug 



2l/T' 



3~3J^S+p — 4Us—2l^T 



r=l I 



r(7 + p-4:/,-2i.T) 



2p 14 D p 7 ' 

' "3 r ^^'^T ~ "3 ~ 3 ^7+p-4i/s -2i/T 



„a6c _ 
"1 — 



/i„ ft,: 



a6c 



n 



abc 



3 

'abc 



2+{^-'^s][-i^ + ^] + [^-K]k-i{h + kc){^ + ^ 



3 



kc_^kb 
kb kc 



,7^5*"^ — ka{kc + Afe)) 



— i's \ — iCs ^ Crp ^ Kd^^" , d^'^ — {kb "I" "I" hkc + kakb + kakc),d'^'^ — ikahkc 



kbkc 



1^8^ {kbkl + kaikl + k^c+ kbkc)} 



,d5 



•abc 



kukf, -\- ka ( kukc ~\~ kbk 



(3.63) 



For both DS and BDS limit the overall normalization factor for three types of polarization can be expressed as: 



1 



-POL 



:u=l(E-mode) 



= { \& :u=2(E<^B mode) 
256 :u=3(B-mode). 



(3.64) 
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As mentioned in the previous sub-section, performing basis transformation cross bispectrum for two scalars and 
one tensor can be expressed as: 



(C(^i)C(fc2)0^fc3)) = {27rfS^^\h +k2 + ks)Bi^'^''\h, k2, ks). 
where we have used the following parameterization: 



B 



ZocoH^;-^ p2 

iVL;3l ^u- 



In this context the polarized non- Gaussian parameter for two scalar and one tensor mode[/; 
as: 



(3.65) 



(3.66) 



NLrz\ ^'^ can be rewritten 



NL;3\x 



where in the DS limit the dimensionful coefficients 



f r 1^ 



1^' 



1^' 



(fc2,fcl,fc3)/C(«)(fc2,fcl,fc3)} 

(3.67) 

(fci, k2, ks) can be re-expressed as 



7a 



X = — oi2 -k2- fc3)(fci + k2- ks){ki -k2 + ks), 



^(2) _ ^(1) ^(3) _ 1 (1) (4) _ fc| (1) (5) _ ,2 (1) (6) _ _§_ (!) 

7a -7a ' 7a - .2 7a > 7a - .2 7a > 7a - %7a > 7a - 1.21.2 7a • 

Most surprisingly, the above coefficients are independent of A due to no parity violation. 
In the EDS limit we have: 



(3.68) 



r fiocan" 



BDS 3Etifc? 



ct-'^-^4'^-\kik2)-''k 



■3 



X 2 



i4!ya+2i/T-12 



r(^^.) 



r(i) 



r(j^T) 



r(i) 



(l-ev-4)'(l-ey-a^)Vi(,/>) 



PL 



where in the BDS limit we have: 

{ h.\ ux" I la' la" \ i,a' i,a" I la' la" \ i,x' i,x" I i,a' jla' 



(3.69) 



c, (i - it) 



O^A,A„, 



' (A:2^' A:3^" P123 + k^' kf P132) (fc^ ^3^" P213 + ^^3^' kf P231) (^i^' kf P312 + k^' kf P; 
(A:^' A:3^" i?i23 + k^^'kf R132) (k^' k{ i?2i3 + k{ k{ R231) (k^' kf R312 + k^ k{ R321 

^ k^nMk2y^ 



fcr(A;2fc3)"= k^^ikiks)'' 

' kf 7?i23 + fc3^" ^2^" ^132) ^ (fcl^' fc^" -R213 + ^3^' kf R231) ^ (k^' k( i?312 + k2k^' R321) 

A;r(A;2/e3)''= ' k^2^ {kik^)-^' ^ ^ ^^(fcife)'^ 



i5a,a„, 

<5a,a„ , , 

<5a/A// 5 



/III II II \ /III I ll\ fill I II \ 

/ I.A i.A I I.A I.A \ / h.X i.A I I.A I.A 1 / i.A i.A _i_ I.A j^A \ 

I fi,2 "-S T n-S n'2 ) I "'I "-S T /C3 I 1 AC]^ K2 T «;2 ""l J 

•2-^ ^ 4- kS-^ i- + k'i- 



fcr(fc2fc3)"= ' k'i^ik^k^)-' 



k^,^{k,k2)-' 



kl 



OSxiXii, 

[k^ Li:2^ + k^' kf Z132) ^2 (fci' l21:i + /m' /-i" i2:il) ^ ^ (^'f /.2" i:ii2 + ^ o' ^'i " ll21 



fcr(^2fc3)^= 



k2^{kiks)-' 



■kl- 



kl''{kxk2)"^ 



where k} = ki where i = 1,2, 3. 



6xi; 



(3.70) 
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In the equilateral limit the expression for the non-Gussian parameter (/jvz,) reduces to the following form: 



r f local 1 
\_JNL;Z\ 



DS 



9fc3 Ytc% ^ 

J q=l 



L/AfL;3j A 



2QLP0L5^^, ((2c, + ct)/c)^'^=+2---9 [Cos ([z. 



{k,k,k)lC^'i\k,k,k) 
)]-[Cos{[ur-\] l)\ 



(3.71) 



1] TtM 3 

2\ 2 



9A;3 

^ / 241^3+21^-12 



r(i) 



r(:^T) 



r(i) 



(l-ey-s^)"(l-6v-4)Vi(0) 



PL 



equil 
X' 

(3.72) 



D. Three tensor correlation 



The interactions involving three tensors are given by the following third order perturbative action: 
('^^''Ofc/jfc ^ / ^^^^^ ^Y^hijfijkhki + ^^^2 {^ikhji - ^hijhki^ hij^ki 



(3.73) 



Now following the prescription of in-in formalism in the interaction picture three point one scalar two tensor correlation 
function both for £'5 and -6-05 can be expressed in the following form: 



(3.74) 



{hi,j,{ki)hi^j^{k2)hi^h{k3)) = -i dr] a (0| \hi,j^{ki)hi^j^{k2)hi^j:,ik3), ([Hint{v)]nj,i,j,i,jX^'\ |0) 

J —oo finri- 

= {2nf5(^\k + fc2 + k) {Bf./.4i,i,i.,.i3i3 (^1' ^2, fc3), 
where the total Hamiltonian in the interaction picture expressed in terms of the third order Lagrangian density as 
{\.Hint{v)]i^j^i^j^i^j^ = - 1 d^^i(^3)hhh\i,j,i^j^i^j3- In this context the bispectrum for three tensor correlation 
can be expressed as: 



{Bhhh}i^j^,^j^,^j^ iki,k2,k3) 



A 

3 ■^«Ul«2j2«3j3 



r rlocal 1 p2 
[JNL;4i 



(3.75) 



where the symbol ; 4 represents three tensor correlation. In this context the non-Gaussian parameter for DS and BDS 
limit can be expressed as: 



r flocal 1 



POL 



A''iJiJ2j2»3i3 (^1' ^2, fcs) + -4f^jti2j2«3j3 (^1' ^3, ^2) + A'^Ji ^2 J2 »3 J3 (^3, ^2, ^1 



+ -^Si»2j2i3j3(^3, ki,k2) + ■A^l^jli^j^i^jS^^, ki,k3) + Ai^jli^j^i^j^{k2, ks, ki) 



where in DS limit we define 

•^niii2j2i3j3 ('''I' ^2) fcs) := 



^ ^ ^ JJcr k k k \ ^ ^ ^ 

'- '""^"^^ ^^3^ 2-/ -^'iJi-'"'(^'l)-^^2j2,m«(fc2)M3j3,»f(A:3) 



l.m.ri 



K 

16 



k3kk3lJ^i3j3,ij{h) - -^k3ik3kMi3j3^jl{k3) 



In the -B-D5 Umit 

r Aocal'\ 



»3J3 



(fciA:2A:3)2-- 



r(^^.) 



r(f) 



r(i/T) 



r(i) 



3 3 9 9 3 



Vp=i 



(3.76) 



1 - E ^'^i - {M.,-„ifc(fci)M2.2,.Kfc2) (3.77) 



(3.78) 
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where K = ki + k2 + ks and the polarization index u = 1{E — mode),2{E^ B — mode),3{B — mode). In BDS 
limit we have: 

Kihhish = Y2^nn-,ij^i2huk^fish-Mcl - ^t^ [M123 + M132 + M213 + M231 + M312 + M321] , 



with 



A 



(2) 



.(3) 



■J^iimik-N'iimji-^hmij [hkksi + k2kk2i + kikku] Q, 

T 

Yt 



»1J1»2J2»3J3 



= -2^-^«ui;»3j3-^«2j2;W [kskksi + k2kk2i + kikku] Q 



(3.79) 



(3.80) 



p=l 



bf'= = (I - jyj,) , bf" = bf^K, 

kl{kb + k^) - i/T^ + {ka{kl + kl + kbkc) + kbkcikb + k^)} Q - i^t 
il - {^li^l + + kbkc) + klkl) + Q - kakbkcikb + k^) 

bf' =(l- Z^t) kakbkc [kbkc + ka{kb + fee)] , hf^ = -klklkl 



bf" 



"5 — 



For both DS and -6.05 limit the overall normalization factor for three types of polarization can be expressed as: 

Wr*^"^ = < 64 :u=2{E^B mode) (3.81) 



4 :u=:l(E-mode) 
,POL _ ) g4 :u=2(E(^B mode) 
1024 :u=3(B-mode). 



After performing basis transformation the relevant three point correlation function for three tensor interaction can 
be expressed in terms of bispectrum as: 



(0''(^i)0''(fc2)0''(fc3)) = {2n)H(^\h+k2 + h)B'i', 



hhh 

A2 , A3 • 



where 



r,hhh _ (27r)'^P^ .Ai,A2,A3 6 



jAi,A2,A3 _ '± fj-iocail" p: 

Ai,A2,A3 - -rrS i,z-^(C,Ch) ~ 5 L-'^i;4JAi,A2,A3 ' 



.2 
u 1 



(3.82) 



(3.83) 



where the the non-linear parameter is given by: 

r fiocan'^ 



3^f(Aifci,A2fc2,A3fc3) 



4Yt 



+ 



K 
32 



kik2k3 



(Aifci + A2A;2 + A3A;3) 



(3.84) 



where in DS limit we define a new function as: 



F{a,p,p) ■.= l-^l^{a + p + p)^{a-p + p){a + p-p){a-p-p). 
o4 a^p^p^ 



(3.85) 
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Similarly in the BDS limit 



BDS 3E,tlfc? 



fl-ev- 



r(i) 



r(i) 



1-ey 



3 3 9 9 3 



where in BDS limit the helicity dependent functions are given by: 

'3 



(3.86) 



A 
A 

A 

with 



(1) 

Ai A2 A3 
(2) 

Ai A2 A3 
(3) 

Ai A2 A3 



Yt 

2^2 AiA A2A A3A 



AAA 



123 ^ -^-^132 



AAA 



M. 



AAA 



A2 A3 



AAA 

312 



M3V ^ 



Ft 



A ;.A 



.A , A 



2c; 



r'^A"A"^A2Ai<5A3A" 



^3 ''■a 



^2 ^2 



(3.87) 



' a'a"a"' r(6+p-6:.T) gA'" ^^QA"'r(9+p-6z.T) 



(^3 



abc' 



p=i 

a' a" a'" / 3 



AAA 



= 1 



Labc 



AAA 



p=l 
^ a'a"a"' 

II 

'A n2 



A \2 , , A , A 



.A \2 



AAA 



{kif{ki )U^^^u{k^{{k^ r+{k^ f+kt k^ )+k^ k^ {kt +k^ )} 



{ (fca' nikf f + ikf f + k{ kf ) + (kf kf )2 } + Q _ fcA' fcA" j^f (fcA" ^ fcA'" ) 



-1,02^ =iK^ ,0^ ^-{k^ kb + k^ kc + k^ ka),0^ =-ik^ hk,, 
where ki — X ki, k2 — X k2 and = A ^3. In the equilateral limit we have: 



A2 ,A3 



Ds 31104(AiA2A3)2 



— + y(Ai+A2 + A3)2 



(3.88) 
(3.89) 



BDS 9k^ 



11 



r(^^.) 



r(|) 



(1 - ey - 3^)'" (1 ~ ey - Sy) V^(0) \ / A . (p);e9m; 



3 9 9 ^ 



f 3 ^ 

^ ^ ^AiA2A3 
Vp=l 



(3.90) 



Numerical values of all such non-Gaussian parameters from three point correlation for different polarizing modes 
are mentioned in the table(|). In this context PC and PV stands for the parity conserving and violating contribution 
for graviton degrees of freedom. 

IV. TREE LEVEL TRISPECTRUM ANALYSIS FROM FOUR SCALAR CORRELATION 

To derive the expression for scalar trispectrum for D3 DBI Galileon let us start from fourth order perturbative action 
up to total derivatives. Consequently the fourth order perturbative action in the uniform gauge can be expressed as: 



(4.1) 
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where the co-efficients t(i{i = 1, 2, 3) defined as: 

_^ Kxx 



as: 




= 




U2 = 


(fKxxx 




kxx 




2 



+ kxx), (4-2) 



Using in-in procedure the four point correlation function both for DS and BDS can be expressed in the foUowing form: 

3 ,0 



(c(fcl)c(fcl)c(fcl)c(^l)) = -iY.I dr^a (oi c(fcl)c(^"^2)c(fcl)c(fc4), (hIi^m),,,, 

j—i J —00 L 



(4.3) 



= (27r)=^^(3)(fc; + fca + kz + k4)Tc{ki,k2, ki), 



where in the interaction picture the total Hamiltonian can be written in terms of the fourth order Lagrangian density 

as: {Hint{n)) ffff = X]?=i ( -'^iiit(^) ) — ~ I d^xCi. Here the trispectrum 7^(A;i, ^2, fcs, ^4) is defined as: 

^^^^ J \ /cccc 



(27r)6p3 



Tc{ki,k2,k3,ki) = —4 -M^ikxMM^'k^) 

11^=1 

= tt4 Ts [(^1^2 + ^I'^'l) (^13^ + + (^1^4 + ^2^3) ('''12'^ + ^13^) 

lli=l«i 

4- ^■^3)^3 1 ju3ju3Wi,-3 I 7^-3\] I Jocai p i 54 JocaJ p3 I 

(4.4) 

where A^(j(A:i, A:2, fcs, ^4) is the shape function for trispectrum and 



-Pc(i) = Pciki3)Pdk3)Pdk4) + Pdkis)Pc{ki)Pc{k2) + Pc(fe4)Pc(fe)^c(^i) + ^c(fe4)Pc(fe)^c(^4) 

+ Pdk23)Pdk2)Pdki) + Pdk23)Pdk4)Pdk3) + Pdki4)Pdk4)Pdk3) + Pdki4)Pdki)Pdk2) 
+ Pdki2)Pdki)Pdk4) + Pc(fci2)Pc(fc2)Pc(A:3) + Pdk34)Pdk3)Pdk2) + Pdk34)Pdk4)Pdki), 



(4.5) 



^C(2) = Pdk2)Pdk3)Pdki) + Pdki)Pdk3)Pdk4) + Pdki)Pdk2)Pdki) + Pdki)Pdk2)Pdk3) 
such that 

+ (k-ik^^ + ^2^3) (A;]^2 + ^13 ) + (^1^3 + ^2^4) (^12 + ^14 )] -'^c(2) 

and t]^£''' and ' ^""^ ^'^^ '^'-'^ linear parameters which carry the signatures of primordial non-Gaussianities of 
the curvature perturbation in trispectrum analysis. By knowing rj^£"' the other parameter 5^^"' can be calculated 
by making use of the following relation: 

9nT^ — [(^1^2 + ^3^4) (^13^ + ^14^) + (^1^4 + ^2^3) (^^12^ + ^13^) + ('^'l^S ^4 7^ 

4- i.3 1.3-1 { 1,-3 I t.-3\] local 

where K = k\ + k2 + k^ + ki. So, there is only one independent piece of information, namely rj^™', that carries 
information about trispectrum. 

To proceed further we denote the angle between ki and kj (with i ^ j) by 6jj then 



Cos(ei2) = Cos(e34) 

Cos(e23) = Cos{Qia) 
Cos{Qi3) = Cos{e2i) 



Cos(03), 

Cos(ei), (4.8) 

: Cos{e2) 



subject to the constraint 

CosiQi) + Cos{e2) + CosiOs) = -I (4.9) 
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comes from the conservation of momentum. 

The exphcit form of r^£''' characterizing the trispectrum can be expressed for our model as: 



[(^1^2 ^3^4) (^13 "I" ^14 ) "I" (^1^4 + ^2^3) (^12 "I" ^13 ) "I" (^1^3 ^2^4) (^12 "I" ^14 )] 



16y4c8 



+ 



+ 



k1k'l{k2-ki) 



kik^ikz-ki) 



3(A;3 + ^4) 12fc3A:4 



^ _^ 3(fc2 + fc4) _^ 12fc2fc4 



+ 



^ ^ 3(fci + ^4) _^ 12fcifc4 



if 



klkl{k2.ki) 
klkl{ks.ki) 



^ ^ 3(fc2 + fcs) ^ 12fc2/C3 



K 



+ 



klkljk^M) 
+ {ki.ki){k2-kz) 



3(A:i + k^) 12fcifc2 
1+ - + ^2 



1 + 



K 

'l2i.j{i<j) = l ^i^j 



+ 



hi 



^ ^ 3(fc3 + fcl) ^ 12fc3fci 



K 



(4.10) 



[(fci.fcl)(fclfc4) + {ki.ks){k2-k4) 



+ 



3kik2k3k4 



1 \ 12kik2k3k4 



l^locall 



K 

2«--VCo,s([..-i] I) 



\BDS 



[(^1^2 + ^3^4) (^13"^ "t" ^14"^) + (fcl^l + ^2^1) (^12'' + ^13'') + (^1^3 + ^2^4) (^12'^ + ^14'^)] 



r(i) 



F4cl2(fcifc2fc3fc4)2 



13 



[r(17 - 8i^,)if®Gi - ir(16 - 8z/,s)if^G2 



r(15 - 8iys)K''Ga - iT{14 - 8vs)K^Ga + r(13 -J>Vs)K'^G^ - iT{l2 - 8v,)K-'G(, 



r(ll - 8i^s)K^G7 - iT{10 - 8iys)KGs + Gg 
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ika.ki)i{'S, 4; 1, 2) + (fe.fc4)X(2,4; 1, 3) 



+ (fc2.fc3)I(2, 3; 1, 4) + (fci.fc4)2:(l, 4; 2, 3) + {ki.k2)I{l, 2; 3, 4) + {ki.ks)I{l, 3; 2, 4) 



+ 



(/Cl./c2)(fc3-fc4) + (fci.fc3)(fc2.fc4) + (fci.fc4)(A:2.A:3) 



Zir(13-8i^s) 



.2:2r(14 - 8;^^) _ .2:3r(15 - 8us) _ Z4r(16 - 8vs) 25^(17 - 8us 



(K) 



13 



15 



16 



17 



G2 = iKGu G3 = Gf (Eti fcf) + Gi (Et(.>,)=i 



\ r^i ruj n-r] 



1 hkjkmj + iG* (j2i,j{ijtj) = l H^j^ ' 



Gi 



^1 (nij,m(i>j>TO)=l ^i^k. 
r4 



n 



2 7.2 7.2 



=1 ^j^m^Ti^ 



{jJi,j,m,n(i<j,m<n,i=^m,j=^n) = l kikjkmkn^ , 



where in BDS Umit 

Gi 
G4 
G5 
Ge 
G7 
Gg 
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and 

m, n) = [K^^/G^^{10 - 6z/«) + K^^/G[{ki + kj + km + fc„)r(ll - 61/^) 

-ff ^^TO^n '^G^kmknikm "I" ^n) kmkn \/ G\ fcjfcj \J G\ (fcj + kj)(kjn + ^n) V'Gi 1 r(12 • 

+ ^ + kn)\fG[ - {ki + kj)kmknVG^] r(13 - 6u,) + kmkn^/GlT{U - 6Us) . 



(4.11) 



(4.12) 



= (ntj,m,n(«>j>r„/„) = l ^'fcj^m^") , Gg = (jltl ^f) ' = 1, 

= i.^, ^3 = (ni,j(i>j) = l ^i^j) ) -^4 = (llij,m(i>j>m) = l kikjkmj , -^5 = (ni=l ^i) 



and in both the limit we have used 



A:i4 = ^23 
A;24 = fcl3 
A;34 = fcl2 



kf + kl + 2kik4Gos{ei), 
k1 + kl + 2k2k4Cos{e2), 
kl + kl + 2k3k4,Gos{e3)- 



6i/,, 



(4.13) 



(4.14) 



22 



Further, using the equilateral configuration (fci = k2 = = = k and K = 4fc) and incorporating the contribution 
from the maximum shape of the trispectrum {Cos{Qi — Cos(02) — Cos{Q^) — and kij{for i < j) ~ -^) the 

non hnear parameter tnl can be recast as: 



equil 



'JGX 



DS 48V3efM|,^L2 
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3GX 



2^^3 
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4 1 - TC^ 



1 - -c 




(4.15) 
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BDS 
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- 16384ir(16 - 8t/s)fc^G2 



(1 - - 



equil 



4096r(15 - 8iy,)k^G'r'^"'^ 



13312/c5 

3 



65536r(17-8j^,)fc«G'i«"' 
' 1024zr(14-8z/,)fc5G'f™' 



256r(13 - 8z/,)fc'^G5«"*' - 64ir(12 - 8;/,)fc3Gg«"' + 16r(ll - 8iys)k'^Gj 



•equil 



4ir(10 - 8iy,)kG\ 



equil 
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^equil ,2 



fc^r(14- 8l/s 



268435456 



g: 



equil 
9 



1024A: 



_'j-equil 



838861^./, 



Z|«"'fcr(15-8:/s) 



1073741824 



^ equil 



12 \ 

r(i6 - 8vs 



4294967296 



fc3r(13-8i 



67108864 



Zf™'r(17 



8vs 



17179869184 



(4.16) 
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FIG. 2: |(a)[Variation of the tml vs sound speed for scalar modes cg, (b) Variation of the tml vs tensor to scalar ratio r in DS 
and BDS limit. 



In this context we have used X'^''™' = 1{hj\ in, n)|i=j=m=„(Vi) and the similar argument holds for C^''"*' and z^'^^^^ 
\/i. The numerical value of t^'^*' equilateral limit is obtained from our set up as 38 < t^'^*' < 89 in BDS limit 



and 48 < r^'^'' < 97 in DS limit within the window for tensor-to-scalar ratio 0.213 < r < 0.250 which is significantly 
large from other class of DBI models. Consequently this model can directly be confronted with the forthcomming 
observational data from PLANCK |^. The graphical behavior of t^'"' with respect to the sound speed (cs) and 
tensor to scalar ratio (r) are plotted in the figure(H)(a) and figure(||)(b) respectively. It is obvious from the figure(||) 
is not much sensitive to the tensor to scalar ratio (r) or the sound speed (eg) within a specified range 



that 



equil 



applicable for our model in DS and BDS limit analysis. In the next section, will show explicitly how this sensitivity 
problem is resolved in DBI Galileon. 
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VIOLATION OF MALDACENA THEOREM AND SUYAMA-YAMAGUCHI RELATION 



To show explicitly the violation of Maldacena theorem, we concentrate on squeezed limit configuration (fcj ~ fcm 
kn with i, m, n = 1, 2, 3 for i 7^ m ^ n) in which the non-Gaussian characteristics can be expressed as: 



NL\DS 



24c| 



(1 - nc) -ris - 26gx - 3sy + 



veW) 



(5.1) 



r(2- ^ 



r(|) 
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Ys6v 
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2Fs 

1^ / 4ey-ys(3-ev) 
4c2 1^ 4(1 + ey) 

^Y[l + Ysf{Ys-cv) 

— + 



"-4 



J-3l|'(nc-l) + §ir(^>) 



X|.(^) + 7^ [J-eI«'(n, - 1) 



5 . -^^ 

(i + yv)(y:v-rv j(i-(i 



2ysc2(l + ey)' 



4c2(l + ey)2 



Consequently the three point correlation function in the squeezed limit can be expressed as: 

(C(^0C(fcl)C(fc3)>^^ = l{2^f5^^\kl + k2 + k3)r^\P{ki)P{kn) 

which, in £'5' limit, turns out to be 

1 " 
44 



(5.2) 
(5.3) 
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and in limit it looks 
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Aey ~ yg(3 - ey) 
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ir(z>) + 



J-3l3««(nc-l) + -|X|'(^) 



s 



+ ^xr(^) 



{l-ev-s^y)\l-VYsf{Ys-ev) 



2Yscl{l + ev)3 

Additionally wc have used in BDS limit the following expressions 
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fel 
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[r(l + x)+r(2 + a;)(i + (3 + a;)^ 



r(i + x) + r(2 + .T)^]}, 

{^^ + (fc;.fc;)A:i}, 
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Here 



(C(fcOC(fcl)C(fc3))r/l^^ = (2^)3<5(3)(fc; + k2 + fc;)(l - nf^)PDs{h)PDs[kn) 



(5.7) 



(c(fc;)c(fc;)c(fci))i^ 



= (27r)3<5(3)(fc; + /Ta + fc3)(l - nf''^)PBDs{h)PBDs{kn) 



(5., 



are collectively the well known Maldacena theorem which is widely accepted almost as a no-go theorem in calculating 
non-Gaussianity. However our result shows significant dev iation from Maldacena theorem. Thus the most appealing 
and strong feature of our model coming from equation ( |5.3| ) is that the violation of well known Maldacena theorem at 
the tree level of non-Gaussianity. This essentially means one can have large non-Gaussianity even from a single field 
model of DBI Galileon. This feature is demonstrated by the non-linearities in the plots in figure (^. This directly 
confirms the credentials of DBI Galileon as the old DBI models fail to give rise to this result. 
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FIG. 3: Variation of ln{fNL) vs ln{TNL) in the equilateral limit configuration in DS and BDS limit. 
Next using equation(^?9|), equation(3.£), equation ( 1. 10| ) and equation(4.11) in DS and BDS limit we get: 
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Jf := lim . 
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where we use the low energy IR momentum cut-off in computing the integrals. This leads to 

IF' > I?'. 

jBDS > jBDS^ 

Substituting the explicit form of the bispectrum and trispectrum in equation ( ^.13| ) we get 



lim 



k\dki 



kidk2Pl'Hki)Prik2) \ [tnl]ds - ^ i[fNL]osy } > 



36 



25 



(5.12) 



(5.13) 



[rNL]DS > ^ iilNLjosf 
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kfdk, / fc2^dfc2Pc (fcl)^C (^2) [rNL]BDS - ^ (\fNL]^r..Y 1> > 

/fcl=fcf« Jfc2=fc^« I 

(5.15) 

This relation directly confirms the violation of standard Suyama- Yamaguchi relation in the context of single field 
D3 DBI Galilean model in SUGRA background thereby solving the generic sensitivity problem of DBI infiation. 
The above results clearly show that the so-called sensitivity problem is being resolved due to the violation of the 
Maldacena Theorem as well as the Suyama- Yamaguchi relation. The most probable reason for such violations are the 
appearance of non-canonical kinetic terms, its derivatives and non-standard stringy frame functions of DBI Galileon 
and Gauss-Bonnet correction in the SUGRA AdSs bulk geometry. Such contributions directly affect the co- 

efficients sitting in front of each an every term in the third and fourth order perturbative action. This results in the 
violation of standard results in non-Gaussianity which can solve the well known sensitivity problem in the context of 
DBI infiation. Very recently other novel aspects of the violation of well known consistency relations in the context of 
single field infiation has been studied in . 

In figure(^) we have shown the parameter space for the non-Gaussian parameters in the logarithmic scale for 
equilateral limit configuration. More technical details can be found in the Appendix. 



VI. SUMMARY AND OUTLOOK 



In this article we have elaborately studied the primordial non-Gaussian features from our proposed model of DBI 
Galileon infiation in D3 brane. We have derived the expressions for three and four point correlation functions in terms 
of the non- linear parameters /nl and tjvl for local and equilateral type of non-Gaussian configurations in the nontrivial 
polarization modes. Hence we have explicitly demonstrated the violation of the well known Maldacena theorem and 
Suyama- Yamaguchi relation leading to a better observational bound for the above mentioned non-linear parameters 
as estimated from WMAP7 [Q dataset. Nevertheless, the results would fit better with the observational predictions 
of PLANCK as expected to be available from their forthcoming datasets. The most significant outcome of our 
model is that it solves the so called sensitivity problem of DBI inflation by introducing Galileon degrees of freedom 
in effective D3 brane in SUGRA background within a specified range of parameter space. This directly confirms the 
physical importance of single field DBI Galileon over the other class of old DBI models as they fail to give rise to such 
promising result. These detectable features lead to the conclusion that this type of models can be directly confronted 
with upcoming data of PLANCK in near future in order to test the credentials of DBI Galileon. 

The remaining open issues in the context of non-Gaussianity for DBI Galileon are studies of mass spectrum of 
primordial black hole formation [|^, as a tool for constraining non-Gaussianity at small scales, effect of the 
presence of one loop and two loop radiative corrections in the presence of all possible scalar and tensor mode fluctu- 
ations in the bispectrum and trispectrum, study of different shapes in equilateral, local, orthogonal, squeezed limit 
configuration for the tree, one and two loop level of non-Gaussianity and calculation of other higher order n-point 
correlation functions to find out the proper consistency relation between all higher order non-Gaussian parameters 
and the analysis of CMB bispectrum and trispectrum in the presence of Galileon in SUGRA background. Given the 
promise the results of the present paper shows, these open issues worth exploring in future as they may give rise to 
even more surprising results. 
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TABLE I: Different non-Gaussian {[fNL;A] ^ ^ ) parameters related to the primordial bispectrum for A=l (three scalar), 2(one scalar 
and two tensor), 3(two scalar and one tensor), 4(three tensor) with polarization index u = 1{E — mode), 2{E^ B — mode), 3{B — 
mode) including all hclicity degrees of freedom represented by Ai, A2 and A3 for DS and BDS limit estimated from our model. In this 
context stands for two projections of helicity for graviton degrees of freedom and "0" represents helicity for scalar mode. 

Here PC and PV stands for the parity conserving and violating contributions appearing in the tree level primordial bispectrum 
analysis. 
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Appendix 



Here throughout the computation we have used the following in-in oscillator modal identities 
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FIG. 4: Tr ee l evel Feynman Diagrams: fa) Th ree scalar correlator, [(b)] Three tensor correlator, (c)One scalar two tensor 
correlator , (d)One tensor two scalar correlator, (e)Four scalar correlator. In this context we have used spiral line for tensor 
contribution and dashed line for scalar contribution. 



The tree level Feynman diagrams for three and four point correlation is drawn in figure (^. 
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